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Abstract 

The observation of the Egret experiment of an excess of diffuse gamma rays with ener- 
gies above E~ = 1 GeV has previously been interpreted in the context of the minimal 
supergravity model (mSUGRA) as coming from neutralino annihilation into mainly b- 
quarks in the galactic halo, with neutralino mass in the vicinity of 50-70 GeV. We observe 
that in order to obtain the correct relic abundance of neutralinos in accord with WMAP 
measurements, the corresponding neutralino-proton direct detection (DD) rates should 
be in excess of recent limits from the Xenon- 10 collaboration. While it does not appear 
possible to satisfy the Egret, WMAP and Xenon-10 constraints simultaneously within the 
mSUGRA model, we find that it is easily possible in models with non-universal Higgs soft 
masses (NUHM). In either case, gluino pair production from m g ~ 400 — 500 GeV should 
occur at large rates at the CERN LHC, and a gluino pair production signal should be 
visible with just 0.1 fb _1 of integrated luminosity. The NUHM interpretation predicts a 
rather light spectrum of heavy Higgs bosons with ~ 140 — 200 GeV over the whole 
parameter space which would interpret Egret data. Spin-independent DD rates are pre- 
dicted to be just above 10 -8 pb, within range of the next round of direct dark matter 
detection experiments. 
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1 Introduction 



An abundance of astrophysical evidence points to the conclusion that the bulk of the matter 
in the universe is composed not of Standard Model (SM) particles, but of some unknown 
non-relativistic elementary particle known as cold dark matter (CDM) [I]. An analysis of the 
three-year WMAP and galaxy survey data sets [2] implies that the ratio of cold dark matter 
density to critical density, 

n CDM h 2 = P cdm/Pc = O.llltHis (2^). (1) 

where h = 0.74 ± 0.03 is the scaled Hubble constant. While the density of CDM is becoming 
precisely known, the identity of the CDM particle (or particles) is still a complete mystery. 
Although numerous candidate CDM particles populate the theoretical literature, the WIMPs 
(weakly interacting massive particles) stand out in that their thermal abundance can be cal- 
culated, and is found to be in rough accord with Eq. (CQ) provided the WIMP mass is of order 
100-1000 GeV. Of the numerous WIMP candidates in the literature, the lightest neutralino of 
supersymmetric (SUSY) theories is especially popular because SUSY solves a host of theoretical 
problems associated with the SM, and also receives some (albeit indirect) support from data (in 
the form of the measured gauge couplings unifying at Q = Mqut under MSSM RG evolution 
and also from other precision electroweak measurements [3]). 

There is at present a multi-pronged effort aimed at identifying WIMP dark matter particles 
and measuring their properties [I]. The most direct approach is to try to detect relic WIMPs 
left over from the Big Bang by observing WIMP-nucleon collisions in experiments located deep 
underground. Limits from the CDMS[5] experiment and more recently from the Xenon-10[6] 
experiment have begun probing the upper limits of SUSY model parameter space. 

WIMP particles can also be searched for at collider experiments such as those at the CERN 
LHC, especially if the dark matter particle is but one of a whole family of particles, some of 
which can be produced via strong and electromagnetic interactions. The dark matter particle 
would then be produced by cascade decays of heavier particles, and would lead to missing 
transverse energy in collider events. Such is the case of theories such as .R-parity conserving 
supersymmetry[7], .fCfT-parity conserving universal extra dimensions (UED)[S] and little Higgs 
models with T-parity[H|. 

Dark matter may also be searched for indirectly. For instance, the sun can sweep up WIMP 
particles as it traverses its galactic orbit, so that WIMPs accumulate at a high density in the 
solar core. Then WIMP-WIMP annihilation to SM particles can occur at high rates in the 
solar core. While most SM particles would be absorbed by the surrounding solar medium, 
multi-GeV scale z/^s would escape and later convert to muons in neutrino telescopes such as 
Amanda/IceCube or Antares. 

In addition, dark matter in the galactic halo can be searched for indirectly via its annihi- 
lation into high energy gamma rays or anti-matter. In the case of gamma rays, searches look 
either for WIMP- WIMP annihilation directly to 77 pairs (loop-suppressed since WIMPs are 
electrically neutral) or via WIMP- WIMP annihilation to qq pairs, followed by q — > 7r° — ► 77 
via hadronization and decay. In the latter case, one expects a diffuse spectrum of gamma rays 
with energies E 1 < Mwimp atop a background arising from cosmic ray spallation onto nuclei, 
inverse Compton scattering and bremsstrahlung. 
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The spectral shape of the gamma ray sky has been measured in the 1990s by the Egret 
experiment in the energy range of 0.1-10 GeV[TU], where already an excess of signal with 
E 1 > 1 GeV above expected background was noted. The E 7 > 1 GeV excess apparently is 
seen in all sky directions. An analysis by de Boer et al. [11] explains the Egret GeV anomaly 
as coming from neutralino annihilation in the galactic halo into mainly bb pairs. A fit of the 
neutralino hypothesis to the Egret data favor a neutralino in the mass range ~ 50 — 70 
GeV. On the astrophysics side of the de Boer et al. interpretation, the strength of the signal 
depends on the dark matter density distribution throughout the galaxy. In order to explain 
the Egret GeV anomaly, de Boer et al. invoke a DM density distribution involving two rings 
of dark matter at 4 and 13 kpc|12j. 

De Boer et al. further interpret the apparent WIMP annihilation signal in the context of 
the minimal supergravity (mSUGRA) model p3|, which allows a complete determination of all 
super-particle and Higgs boson masses and mixings in terms of just a few parameters |14] 

mo, m±/2, A , tan/? and sign(fi), (2) 

where m is a common scalar mass at energy scale Q = Mgut, m i/2 is the common gaugino 
mass at Mgut, Aq the common trilinear GUT-scale soft breaking mass, tan (3 is the ratio of 
higgs vevs and /i is the superpotential Higgs mass parameter. The magnitude- but not the sign- 
of the /i term is determined by constraints arising from requiring an appropriate breakdown 
of electroweak symmetry in the mSUGRA model. In mSUGRA, once the GUT scale soft 
SUSY breaking terms are stipulated at Mgut, their weak scale values can be calculated via 
renormalization group equations (RGEs). The physical SUSY particle masses and mixings 
can then be calculated in terms of the weak scale soft SUSY breaking terms via well-known 
algorithms [7]. 

To generate a neutralino WIMP with mass 50-70 GeV, de Boer et al. require m\/2 in the 
range of 130-170 GeV. This rather low value of m 1( / 2 typically leads to light Higgs masses 
rrih ~ 114.4 GeV (the limit from LEP2 on SM-like Higgs scalars) and to large contributions 
to the branching fraction BF(b — > 57). To avoid these constraints, de Boer et al. adopt a 
large value of niQ ~ 1500 GeV[15j. Then, to gain accord with the measured relic density (Eq. 
[Q, they require large tan/? ~ 54. For such a large value of tan/?, the 6-quark and r-lepton 
Yukawa couplings become very large, the pseudoscalar Higgs mass m A falls, and its width grows. 
Neutralinos Z\ can then annihilate efficiently in the early universe as Z\Z\ — > bb via virtual, 
non-resonant A exchange in the s-channel|16j. Just as neutralinos can annihilate efficiently via 
A* in the early universe, so can they annihilate efficiently via A* to bb in the galactic halo, since 
the A-annihilation diagrams are s-wave (while h, H annihilation is p-wave, and thus velocity- 
suppressed). The dominant halo annihilation Z\Z\ — > A* — > bb — > ir° — > 7 can then describe 
the Egret gamma ray excess. Note that there is little uncertaintly in the gamma spectrum 
from 50-70 GeV WIMP annihilation, since the corresponding process e + e~ — > bb — > 7 has been 
well-measured at LEP/LEP2. 

It is important to note that several alternative explanations/insights regarding the Egret 
excess have emerged [T7]. 

• In Ref . [18] , Strong, Moskalenko and Reimer verify that a "conventional model" of cosmic 
ray production and propagation is insufficient to explain the Egret GeV anomaly, even 
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if augmented by hard sources of additional cosmic ray production in the inner galaxy. 
However, by suitably adjusting the spectrum of cosmic protons and electrons, they find 
they are able to by-and-large match the Egret data. 

• Alternatively, Stecker et al.^9\ claim the Egret excess can be explained by a calibration 
error in Egret gamma measurements with E 1 > 1 GeV. This would explain why an excess 
of high energy gammas comes from all sky regions. 

• Bergstrom et al. [2D] point out that if they use the deBoer derived distribution of galactic 
dark matter, including the ring structure, then the SUSY region favored in the deBoer 
analysis yields an anti-proton flux far in excess of measurements from the BESS experi- 
ment. These calculations adopt the isotropic DarkSUSY model of cosmic ray propagation. 
DeBoer et al. counter that using an anisotropic model of galactic cosmic ray propagation 
would greatly reduce the expected p flux[2T]. 

In this note, we wish to examine the SUSY interpretation of the Egret excess, and compare 
it with other constraints on sparticle masses. We find that the SUSY interpretation in terms of 
the mSUGRA model is in conflict with recent results on direct detection of dark matter from 
the Xenon-10 experiment. By moving to models with non-universality in the Higgs sector, 
however, one may preserve the SUSY interpretation of the Egret gamma ray excess, while 
staying below bounds on direct detection of DM. In any case, the imminent turn-on of the LHC 
should decide the issue. The deBoer interpretation predicts gluinos with mass rrig ~ 400 — 500 
GeV. Production cross sections for gluino pair production in this mass range at the LHC are 
at the 10 5 fb level. Thus, LHC should decisively test the deBoer scenario with as little as 0.1 
fb _1 of integrated luminosity. If the mSUGRA interpretation is correct, then the heavy Higgs 
boson A should have mass around 200 — 400 GeV. If an interpretation in terms of the NUHM 
model is correct, then should be much lower- in the 140 — 200 GeV range. 

2 Confronting SUSY interpretation of Egret GeV anomaly 
with Xenon-10 direct dark matter search 

2.1 mSUGRA analysis 

We begin by calculating sparticle mass spectra using the Isajet 7.76[22] Isasugra code. Isajet 
begins with weak scale DR values for the three gauge couplings and Yukawa couplings, and 
evolves up in energy to determine Mqut, defined as the Q value where gauge couplings g± = gi- 
At Mqut, soft SUSY breaking boundary conditions are input, and the set of 26 coupled 2- 
loop RGEs are evolved down to M wea k- Beta-function threshold effects are included in the 
1-loop portion of RGEs for gauge and Yukawa couplings, giving a smooth transition between 
MSSM and SM effective theories. All soft terms which mix are frozen out at scale Msusy — 
y/rrii L mi R , while all non-mixing soft terms are frozen out at their own mass scale {e.g. rn\ R 
stops running when Q = rriu R is reached[23]. The RG-improved 1-loop MSSM scalar potential 
is minimized at Q = Msusy, which determines the value of /i 2 . All tree level sparticle masses 
are computed. Once these are known, then all 1-loop sparticle masses are computed, including 
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Figure 1: Plot of reach of various present and future dark matter direct detection experiments compared 
to theory prediction in the mSUGRA model. 

SUSY threshold corrections to m t , nib and m T . The threshold effects alter the trajectories of 
the running couplings, so that an up-down iterative approach is used to calculate all 1-loop 
corrected sparticle masses; the iterations terminate when a convergence criterion is satisfied. 
At this point, all sparticle decay branching fractions are calculated, along with neutralino relic 
density Q^h 2 , a M = (g - 2) M /2, BF(b -> sj), a SI (Z x p) , BF(B S -> fJ, + fJ,~) and (av)\ v ^ , via 
the Isatools package [24"] . The latter quantity, the neutralino annihilation cross section times 
relative velocity, in the limit as v — > 0, is the crucial particle physics quantity needed to evaluate 
various halo annihilation processes. 

Regarding the neutralino direct dark matter detection cross section, we note here that a 
new limit on the spin-independent neutralino-nucleon scattering cross section, asi{Zip), has 
appeared from the Xenon-10 collaboration [6]. This new limit, displayed as the solid red curve 
in Fig. Il[25j . excludes asi(Zip) ~6 x 1CT 8 pb for ~ 60 GeV. We also show in this figure 
the projected reach of several future direct detection experiments (dashed curves), along with 
theoretical predictions from a scan over mSUGRA model parameter space[26] (pink region). 
The reach contours assume a standard local DM density of pdm = 0.3 GeV/cm 3 , and a standard 
DM velocity profile. 

We first attempt to verify the de Boer et al. suggested mSUGRA interpretation using the 
Isajet code. The requirement that ~ 60 GeV means that- in models with gaugino mass 
unification and \i ^> M 2 - the lightest neutralino should be dominantly bino-like. The trick is 
to get a low relic density in accord with WMAP, while maintaining a large value of (cry) |„^o, 
so there is sufficient neutralino annihilation in the galactic halo. The LEP2 constraints that 
i). > 103.5 GeV and ii). mf ~ 95 GeV means that co-annihilation cannot be used to 

reduce the relic density to WMAP allowed levels. De Boer et al. suggest taking m^j ~ 160 
GeV to get ~ 60 GeV, and large m to suppress SUSY contributions to BF(b — > 57) 
and to raise the value of to LEP2-allowed values^- These input parameters, along with 

^^Here, we require m/j ~ 111 GeV to account for a roughly 3 GeV slop in the theory calculation of m/ l; while 
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Figure 2: Plot of a). VL~ h? versus tan/? and b). asi{Z\p) versus tan/3 for two cases in the mSUGRA 
model. 

A = and tan/3 = 10, give ^^Ji? ~ 10, which is two orders of magnitude higher than Eq. 
flTJ). In Fig. [2]a., we plot the value of Q^Jt? versus tan/3 for m = 1500 GeV, = 160 

GeV, A = and m t = 175 (as in de Boer et al. [13]), arid m t = 171 GeV (the central value of 
m t as recently measured by DO and CDF[2H3)- While ^^Ji? is too large for most of the range 
of tan /3, we see that it drops to the measured value around tan (3 = 52 — 55. At this high a 
value of tan /3, the b and r Yukawa couplings become very large, while the value of drops 
(see Fig. [3]). Even though 2m^ is still far from the A resonance, annihilation through the 
virtual A* becomes dominant enough to lower the relic density to WMAP- allowed values. The 
beauty of this approach is that when neutralino annihilation in the early universe via s-channel 
A exchange is large, so also is halo annihilation of neutralinos[29]. This contrasts with the case 
of Z\Zi — > h —* bb, where early universe annihilation can be large, but (av)\ v ^o — > 0, so that 
the neutralino halo annihilation rate is small. 

In Fig. we show the spin-independent neutralino-proton scattering cross section versus 
tan/3 for the same parameters as in Fig. [2]a. We see that as tan/3 increases, the value of 
<Jsi{Zip) also increases. This is due to Higgs exchange direct detection scattering diagrams, 
and the increasing magnitude of the 6-quark Yukawa coupling. The important point to notice, 

LEP2 requires a SM-like h to have m h > 114.4 GeV [2?]. 
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Figure 3: Plot of sparticle masses versus tan/3 for the two cases of mSUGRA model displayed in 
Tabled 

however, is that the value of asi(Zip) has increased into the Xenon-10 excluded region at a 
tan j3 value somewhat lower than that needed for the relic density to enter the WMAP-allowed 
dark matter band. Thus, technically, these cases for a SUSY interpretation of the Egret GeV 
anomaly would be excluded by the new Xenon-10 limits. 

We list two WMAP-allowed mSUGRA points in Tabled] The first, suggested in the de Boer 
et al. analysis[l3] has m t = 175 GeV so that an exact comparison can be made to Ref. [T3j . 
The second point has m t dialed down to 171 GeV in accord with recent top mass measurements 
at the Tevatron. In addition, for the second point, we take A = —900 GeV which raises the 
value of rrih more closely into accord with LEP2 Higgs search limits. We see that both cases 
have nig ~ 60 GeV, while m- g ~ 470 GeV. Squarks and sleptons have masses above the TeV 
range. We also list in the Table the relic density Q^^k 2 , BF(b — > 57), the SUSY contribution 
to the muon anomalous magnetic moment Aa M , the branching fraction for B s — > fi + fi~ decay, 
asi(Zip) in pb, and (av)\ v ^o- The neutralino-proton spin independent scattering cross section 
is ~ 3 x 10 -7 pb in both cases- well above the Xenon-10 limit. The value of (av) \ v ^o ~ 2 x 10~ 26 
cm 3 /sec, a value which gives sufficient halo annihilation in models assuming the de Boer DM 
halo profile. 

We next wish to check if Xenon-10 would exclude all mSUGRA interpretations of the Egret 
data. We do so by scanning over the entire mSUGRA model parameter space: 

100 GeV < m < 4000 GeV 

10 GeV < m 1/2 < 1000 GeV 

-3000 GeV < A < 3000 GeV 

1.1 < tan/3 < 60 

fj. > 0. (3) 

Our results of this scan are plotted in Fig. HI Here, we keep only points with 0.09 < fi^/i 2 < 
0.13, and also 50 GeV < < 70 GeV. Green points have a BF(b — > 57) in close accord 
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1022.1 
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m- eL 


1501.6 




1501.6 


853.6 


mg R 


1499.8 




1499.8 


802.4 




110.9 




106.3 


131.7 


m z 2 


110.4 




106.7 


131.0 


^1 


62.4 




61.8 


66.6 


m A 


309.1 




347.0 


157.0 


m h 


113.6 




112.8 


116.6 


BF{b sj) 


0.11 




0.11 


0.10 


3.0 x 10' 


-4 


2.4 x 10" 4 


3.1 x 10~ 4 


BF{B S -> 


10.3 x 10~ 


-10 


10.0 x 10~ 10 


5.4 x 10" 10 


2.2 x 10" 


-9 


9.3 x 10~ 9 


3.7 x 10- 8 


cr sc (Zip) [pb] 


3.2 x 1(T 


-7 


3.1 x 10~ 7 


2.6 x 10^ 8 


(a^)|^o (cm 3 / sec ) 


2.0 x 10- 


26 


2.3 x 10" 26 


1.6 x 10" 26 



Table 1: Masses and parameters in GeV units for three Egret-motivated benchmark points 
using Isajet 7.76. 
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Figure 4: Scan of mSUGRA parameter space for models which obey LEP2 constraints and have 
0.09 < £1^ h 2 < 0.13. We plot results in the (av) \ v ^o vs - °~si(Zip) plane. Green dots have good 
BF(b — > sj) while red dots deviate from the measured branching fraction. 

with the measured value: BF(b — > sj) = (3.55 ± 0.26) x 10 -4 from a combined analysis [30] 
of the CLEO, Belle and BABAR experiments. Yellow and especially red points give branching 
fractions further from the experimental central value. The surviving points are plotted in the 
{av)\ v _>o vs. asi(Zip) plane. We see that most points populate the low (av)\ v ^o and low 
crsi(Zip) regions. These points come from either the stau co-annihilation region or the h- 
resonance annihilation region of the mSUGRA model. There are no hyperbolic branch/focus 
point (HB/FP) contributions since < Myy, and so ZiZi — > W + W~ (as is enhanced in the 
HB/FP region) is kinematically forbidden. The points at large (o~v)\ v ^,q ~ 10~ 26 cm 2 /sec (so 
that they have a high halo annihilation rate into gammas) also are all above the Xenon- 10 dark 
matter limitl For this reason, it seems an interpretation of the Egret GeV anomaly in terms 
of neutralino annihilation in the mSUGRA model is ruled out. Of course, one way out is to 
assume we live in a local void of dark matter, and the local density is far below the assumed 
value of 0.3 GeV/cm 3 . One may also assume much lower local WIMP velocities, which would 
also lead to lower detection rates. (It may also be the case that we live in a locally overdense 
region, or that the velocity profile is harder than expected, leading to larger than expected DM 
detection rates.) Here, we will not further entertain these possibilities. 

2.2 NUHM2 analysis 

While the mSUGRA model does not seem adequate to explain the Egret GeV anomaly in the 
face of the new Xenon-10 limit, other less restrictive supersymmetric models may do the job. 
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One highly motivated model beyond mSUGRA consists of models with non-universal soft SUSY 
breaking Higgs masses (NUHM). In simple SO (10) SUSYGUT theories, the Higgs supermulti- 
plets live in a 10 of 5*0(10), while the matter supermultiplets live in the 16 dimensional spinor 
representation. Thus, one might naturally expect Higgs SSB terms to have different GUT scale 
masses than matter SSB terms (this is the one-parameter NUHM model, or NUHM1|3T]). In 
SU(5) SUSYGUT models, the doublet H u lives in a 5, while the doublet lives in a 5. In 
this case, both m? Hu and m? Hd can be taken as independent parameters, whereas the matter 
scalars remain unified to m . This is the two-parameter NUHM model, or NUHM2[32J. In 
NUHM2, the GUT scale parameters m 2 H and m 2 H can be traded for independent weak scale 
parameters \x and m A (whereas in mSUGRA, these quantities are derived from the GUT scale 
inputs, mainly mo). Here, we will examine the NUHM2 model, with parameter space given by 

m , mi/2, A , /I, m A , tan/3, (4) 

where we take m t = 171 GeV as usual. Our goal will be to lower tan j3, so that we will diminish 
the direct detection cross section to levels below the Xenon-10 limit. Meanwhile, we wish to 
maintain a large (av)\ v ^o so that we maintain a high rate of neutralino halo annihilations. 
This can be done by lowering m A so that we move nearer (but not directly on) A-resonance 
annihilation. 

Here, we scan over the NUHM2 parameter space: 



100 GeV < 


m 


< 2000 GeV 


50 GeV < 


mi/2 


< 300 GeV 


-3000 GeV < 


A 


< 3000 GeV 


1.1 < 


tan j3 


< 60 


50 GeV < 


/' 


< 1000 GeV 


70 GeV < 


m A 


< 500 GeV, 



while again plotting points which satisfy the WMAP relic density bound and have 50 GeV < 
< 70 GeV. The points are again plotted in the (av)\ v ->o vs. asi(Z\p) plane in Fig. [51 
In this case, we do find a collection of points with simultaneously the correct relic abundance 
and neutralino mass, a high rate of halo annihilation (since (av)\ v ->o ~ 10 _26 cm 3 /sec), and 
with crgi(Zip) below the Xenon-10 limit! While the collection of points does satisfy the Xenon- 
10 limit, note that they do not extend to arbitrarily small values of agi(Zip). In fact, one 
prediction is that if the NUHM2 SUGRA model is to explain the Egret GeV anomaly, then 
agi(Zip) ~ 10~ 8 pb, which is well within range of a number of upcoming direct detection 
experiments, including Lux, Xenon-100, WARP-140 and mini-CLEAN. 

As an example of a NUHM2 point which satisfies all constraints, we list the third point in 
Tabled] a point with mo = 831.8 GeV and m±/2 = 161.2 GeV. Like the mSUGRA models, it 
has a light gluino, with nig ~ 450 GeV. It also has m A ~ 157 GeV, which is generically below 
the mSUGRA models, which predict m A ~ 200 — 400 GeV. The tan/5 value is just 17.6, so the 
b- Yukawa coupling is not so large, and the A-width is rather narrow. 

To show explicitly the range of m§ and m A expected in our scans, we plot in Fig. [6] points 
from a), the mSUGRA scan and b). the NUHM2 scan in the m A vs. m g plane. In both cases, 
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log 10 a xp (pb) 

Figure 5: Scan of NUHM2 parameter space for models which obey LEP2 constraints and have 0.09 < 
Vl~ h? < 0.13. We plot results in the (o~v)\ v ->o vs. asi{Z\p) plane. Green dots have good BF{b — ► 57) 
while red dots deviate from the measured branching fraction. 

surviving points have m~ : 50 - 70 GeV and 0.09 < Q^h 2 < 0.13, BF(B S -> < 10~ 7 

and further, we require (av)\ v ^o > 10~ 26 cm 3 /sec so that there is sufficient halo annihilation to 
explain the Egret excess. The points in the mSUGRA case all have m g ~ 400 — 550 GeV, while 
rriA > 200 GeV. In the NUHM2 case, we further require osi{Z\p) < 6 x 10~ 8 pb. Here, we 
see a similar range of m g is allowed. However, in the NUHM2 case, we always have rriA < 200 
GeV. This is an important distinction between the two interpretations which can be directly 
tested/measured at LHC. 

3 Egret SUSY interpretation and the LHC 

The CERN LHC is expected to turn on in mid-2008, and gradually begin accumulating data 
into 2009. It is not unreasonable to expect of order 0.1 fb _1 of integrated luminosity in the first 
full year of running. 

If the SUSY interpretation of the Egret GeV anomaly is correct, then we expect the DM 
particle to be a neutralino of mass ~ 50 — 70 GeV. Most SUSY models also assume 
gaugino mass unification at the GUT scale. In such models, the gluino mass is expected to be 
m g ~ 7 Mi, and if the Zi is dominantly bino-like, with \i ^> Mi, then we expect ~ Mi and 
thus would predict m g ~ 350 — 500 GeV. Squarks and sleptons can be considerably heavier, as 
is the case in the points listed in Table [TJ In this type of scenario, we would expect new physics 
events from SUSY at the CERN LHC to be dominated by gluino pair production, followed by 
gluino cascade decays[33]. The cross section for pp — > ggX at the LHC for the above range of 
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Figure 6: Scan of a). mSUGRA model and b). NUHM2 model parameter space for points which obey 
LEP2 constraints and have 0.09 < ^^Ji? < 0.13 and (av)\ v ^o > 10~ 26 cnt'/sec. We plot results in 
the rrig vs. m,A plane. Green dots have good BF(b — > S7) while red dots deviate from the measured 
branching fraction. 

gluino masses is ~ 10 4 — 10 5 fb[3l]. Thus, with just 0.1 fb _1 of integrated luminosity, we can 
already expect 10 3 — 10 4 gluino pair events to be recorded per LHC experiment. 

While new physics may be lurking in the LHC already shortly after turn-on, it is unclear 
whether the detectors will be fully calibrated to allow for a new physics search. For instance, 
traditional SUSY searches rely on a E™ 1SS + jets signature, where large E™ 1SS is required to 
reject SM background events from the SUSY signal. To properly use the E™ 1SS variable, a full 
knowledge of the detector is required, since E™ ms can not only arise from signal and background 
events, but also from 1. dead regions of the detector, 2. "hot", or mis-firing calorimeter cells, 
3. cosmic ray events and 4. energy mis-measurement in active calorimeter cells. In Ref. [35], it 
was recently pointed out that early discovery of SUSY at the LHC was possible without using 
_E™ 1SS , and that a reach in m- g to 600 — 700 GeV could be attained without using E™ 1SS , and with 
only 0.1 fb _1 of integrated luminosity. Effectively, the idea was to make use of multi-lepton 
production [36] in the lengthy sparticle cascade decays. Thus, requiring > 4 jets production 
along with > 2 or 3 isolated leptons, allowed for an excellent rejection of SM background 
(dominated by ti production) compared to signal so that a SUSY discovery could be made. 
The values of rrig expected from the SUSY interpretation of the Egret GeV anomaly fall well 
within this "early discovery" range. 

Here, we use the same detector simulation, jet finding algorithm and lepton isolation cri- 
terion as detailed in Ref. [35], and adopt the same set of SM background events from QCD 
jet production, W + jets, Z + jets, ti production and vector boson pair production. We re- 
quire first that signal and background events satisfy the set of cuts CV: i). n(jets) > 4, ii). 
E T (jl, j2, j3, j4) > 100, 50, 50, 50 GeV, respectively (where jets are ordered according to E T 
value) and in), transverse sphericity St > 0.2. We then plot a multiplicity of isolated leptons 
(a lepton i = e or /i is isolated if it has E T (£) > 20 GeV, \rj(£)\ < 2.5 and £ E^ Us < 5 GeV in 
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Figure 7: Distribution in number of isolated leptons from the NUHM2 point in Tabled after cuts Cl' . 
We also show component SM backgrounds and the summer SM background (gray histogram). 



a cone of AR = \/ Ar/ 2 + A<ft 2 < 0.2 about the lepton direction). The results are shown in Fig. 
[3 Here, we see that SM background dominates signal for m — or 1. But already at ne = 2, 
the signal from the NUHM2 point in Table [1] stands out above background. At ri£ = 3, there 
still remains ~ 400 fb of signal, while BG is negligible. Given these results, it seems that the 
SUSY interpretation of the Egret GeV anomaly should be easily testable at the LHC after only 
0.1 fb -1 of integrated luminosity is obtained, and even before the detectors are fully calibrated 
such that they can perform a reliable search for E™ 1SS + jets events. 

Once a SUSY signal is obtained, then the set of likely signal events can be scrutinized 
to try to reconstruct sparticle masses, etc. The starting point is often to first look at the 
opposite sign/same flavor (OS/SF) dilepton invariant mass spectrum m(££)^7\. We plot the 
distribution da/dm(££) in Fig. [8] arising from the second mSUGRA point in Tabled], after cuts 
CY and requiring a pair of OS/SF isolated leptons. In this case, a clear mass edge is seen 
at — nig = 44.9 GeV. There is also a Z peak in both signal and BG (the latter arising 
because Isajet includes W and Z radiation in its parton shower algorithm). 

In the case of the NUHM2 model from Table [TJ, the ti mass is so light that g — > tti 
dominantly, and the Z 2 production via cascade decays is somewhat suppressed. Furthermore, 
the Z 2 — > Zie + e~ branching fraction is suppressed to the 0.8% level; in this case, the suppression 
is due to the presence of relatively light A and H Higgs bosons, which enhance the decay 
zT 2 — > Z\bb to the 45% level, at the expense of first/second generation decay modes. Thus, in 
the m(££) distribution for NUHM2 shown in Fig. (9J we see a continuum distribution instead of 
a distinct mass edge (the mass edge would occur at — = 64.4 GeV in this case). The 
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Figure 8: Distribution in OS/SF dilepton invariant mass from the second mSUGRA point in Table 
[7] after cuts CI' . We also show component SM backgrounds and the sum of SM background (gray 
histogram). 



other crucial observable to test between a mSUGRA or NUHM2 interpretation of the Egret 
excess will come from the LHC measurement of the heavy Higgs boson mass spectrum. The 
Atlas and CMS groups have posted reach plots for MSSM Higgs bosons in the rriA vs. tan/? 
plane [381 [39] . The H and A Higgs bosons should be seeable in the bb, tt and even 
modesgD] in the NUHM2 model for tan/3 ~ 6 since uia ~ 200 GeV, and their masses should 
be measureable. The mSUGRA interpretation requires tan (3 ~ 50 and has relatively light ttia 
as well, and should likewise be visible, but with mass rriA ~ 200 GeV. In Fig. [101 we present 
the mSUGRA a), and NUHM2 b). points in the (tan/3, tua) plane with lightest neutralino 
mass in the range of 50-70 GeV. Points obey LEP2 constraints and have 0.09 < Q^Ji 2 < 0.13, 
(&v) >o > 10 -26 cm 3 /sec and Br(Bs —>■ fi + < 10~ 7 . For NUMH2 parameter space we also 
require o~si{Zip) < 6~ 8 pb to satisfy the Xenon-10 limit. One can clearly see that indeed these 
scenarios suggest a very distinctive A-Higgs boson phenomenology. Recent Tevatron limits [JT] 
for the (tan/3, ttia) plane already exclude ttia < 180 for tan/3 > 50. So, the Tevatron can test 
now a small part of NUMH2 parameter space with large values of tan f3, while parameter space 
with 100 < rriA < 200 and tan j3 < 30 could be tested and completely covered only at the LHC. 
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Figure 9: Distribution in OS/SF dilepton invariant mass from the NUHM2 point in Table [7] after 
cuts CI' . We also show component SM backgrounds and the sum of SM background (gray histogram). 



4 Conclusions 

In this paper, we have examined the SUSY interpretation of the Egret GeV anomaly. The 
SUSY interpretation requires m^r ~ 50 — 70 GeV, and a relic density of ^^Ji 2 ~ 0.1. In order 

to satisfy these criteria, while maintaining a large halo annihilation rate for Z\Z\ — >• A* — > 66, 
one must move to very large tan/3 ~ 52 — 55 in the context of the mSUGRA model. At this 
high tan j3, the predicted spin-independent Zip scattering cross section exceeds recent limits 
from the Xenon- 10 collaboration. 

In order to maintain a SUSY interpretation of the Egret GeV anomaly, we suggest moving to 
SUSY models with a non-universal Higgs sector: the NUHM2 model. In this case, with freedom 
to adjust the value of rriA, one can reduce the value of tan/3 so that the predicted asi(Zip) is 
lowered below Xenon- 10 limits, but maintain a valid relic density and halo annihilation rate by 
lowering rriA to values below 200 GeV. If the SUSY interpretation of the Egret GeV anomaly 
is correct, then we predict direct detection cross sections asi{Z\p) ~ 10^ 8 pb, which should 
be accessible to the next round of DD experiments. Further, the gluino mass should be in the 
range 400 — 550 GeV, and should be seeable in the multi-jet plus multi-lepton channel at LHC 
with just 0.1 fb _1 of integrated luminosity. A dilepton mass edge in these signal events may 
not be apparent since the light Higgs spectrum enhances Z2 decay to third generation fermions 
at the expense of first/second generation leptons. Since mA,H ~ 200 GeV, the A and H should 
be readily visible at LHC via searches for A, H — > 66, t + t~ or yU + /i~, even for relatively low 
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Figure 10: Points from a). mSUGRA and b). NUMH2 interpretation of the Egret GeV anomaly 
in (tan/3, itia) plane with lightest neutralino mass in the range of 50-70 GeV. Points obey LEP2 
constraints and have 0.09 < Vl~ h? < 0.13, (av)\ v —>o > 10 -26 cm? /sec and Br(Bs — ► < 10 -7 . 

For NUMH2 parameter space we also require asi{Z\p) < 6 x 10 -8 pb to satisfy Xenon-10 limit. Green 
dots have good BF(b — > S7) while red dots deviate from the measured branching fraction. 



values of tan/3 ~ 6 — 10. This observation clearly distinguishes between the mSUGRA and 
NUHM2 interpretations of the Egret GeV anomaly. 
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